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Abstract. We introduce a method based on the finite size scaling assumption which 
allows to determine numerically the critical point and critical exponents related to 
observables in an infinite system starting from the knowledge of the observables in 
finite systems. We apply the method to bond percolation in 2 dimensions and com- 
pare the results obtained when the bond probability p or the fragment multiplicity 
m are chosen as the relevant parameter. 
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1. Introduction 

Phase transitions are phenomena which are related to infinite systems. Hence 
the quantities which characterize these phenomena must be related to infinite sys- 
tems. In fields like solid state physics this can be more or less easily realized since 
the systems which are considered can most of the time be assimilated to infinite 
systems. But this is not the case in nuclear physics where nuclei contain at most a 
few hundred constituents. 

Among the most striking features which are observed in the study of nuclear 
fragmentation emerges the fact that there are signs for the existence of a phase 
transition. This phenomenon is commonly under scrutiny nowadays, both experi- 
mentally and theoretically. There are indications for the existence of a first order 
transition but its existence is not clearly established up to now [1-4] . On the other 
hand, charge distributions of fragments obtained in peripheral collisions show uni- 
versal properties [5] which are conspicuously compatible with bond percolation [6] 
and phase space models which describe excited and disordered systems [7,8]. If one 
relies on a percolation-type interpretation of these results, it is tempting to conclude 
that one is confronted with the existence of a second order phase transition whose 
characteristic features are smoothed by finite size effects. 

The fascinating perspective that it might be possible to unravel the existence 
of a phase transition in nuclear matter motivated the attempt to extract more 
precise information about this kind of phenomenon. Such information is given by 
the numerical values of critical exponents which allow to fix the class to which 
the observed phase transition belongs. The most important exponents are v which 
governs the correlation length £, a which characterizes the behaviour of 7712/ m\ 
(nik is the moment of order k of charge or mass distributions), (3 which fixes the 
behaviour of mi close to the critical point, r which gives the slope of the power law 
distribution at the critical point and 7 which governs the strength of the singularity 
of 7712 at this point. These quantities are connected by universal relations [16]. 

The determination of critical exponents raises however two intimately related 
points. The first one concerns the obtention of these quantities from observables 
which characterize finite systems and hence do not show any singularity as they 
do for infinite systems. This is a conceptual problem which, to our opinion, is 
fundamental. The second point which is directly linked to it concerns the direct 
extraction of exponents from finite size model results or experimental data neces- 
sarily obtained from systems of finite size L, especially in the case of nuclei. Any 
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such attempt leads inevitably to a size-dependence of any "exponent" say a, i.e. 
a = a(L). As a consequence, it gets difficult if not impossible to compare these 
quantities to corresponding reference exponents characterizing well-known systems, 
like the liquid-gas, percolation or magnetic systems, which have been worked out 
in the infinite limit and fix the universality class to which the considered system 
belongs. The possibility to extract critical exponents from finite size systems has 
been attempted in refs. [9-11] and was followed by controversial debates [12,13]. 

In the sequel we aim to present and use a new numerical approach which has 
been worked out recently in order to analyze finite size calculations in the 2D and 
3D Ising models [14]. The method relies on the finite size scaling (FSS) assumption 
[15]. 

2. The method 

We sketch here the essential steps which lead to the determination of the critical 
point and critical exponents for infinite systems starting from finite systems of 
different sizes L in spaces of any physical dimension. 

Consider an observable Ox, (At) characterizing a finite system of linear size L 
where At measures the distance of the variable t to its value at the critical point t c , 
At = (t c — t)/t c , At > 0. Define x (L,t) = ^(t)/L, where £x,(t) is the correlation 
length in the system of size L for a fixed value t of the relevant variable. Then the 
FSS assumption states that 

Ox-(At) = O (At) • Qo(x (L, t)) (1) 

where Qo{x (£,£)) is a so called scaling function and O (At) fixes the value of the 
observable in the corresponding infinite system. The remarkable property of Qo(%) 
is the fact that it is universal, which means here that it is independent of t for 
different values of x in the range [0, t c ] of the variable t. 

It comes out that quantities like £x,(At) and other observables Ox, (At) converge 
to a finite value with increasing L. Hence, if Ox, (At) is known for different values of 
L and if its asymptotic value with increasing L, O (At), can be numerically reached 
the scaling function Qo(x (L,t)) can be obtained through (1). Once this quantity 
is known for different values of x = £i,(t)/L, it can be used to get 

O (At) = L (At)/Q (x) ~ At" Q - (t c - t)" Q (2) 

where t c is the critical point and a the corresponding exponent. 
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In practice, once Qo{ x ) is numerically known it can be easily parametrized 
in terms of polynomials in x or some other function of x [14]. The knowledge of 
Qo and O^(At) gives O (At) through (2), the numerical result can be fitted to an 
expression like (t c — t)~ a through which t c and a can be determined. The accuracy 
of the fits can be controlled by means of x 2 tests. 

3. Application to 2D percolation 

We have checked the applicability of the method presented above on the 2D 
bond percolation problem on a square lattice for which the numerical values of the 
critical exponents are well known [16]. In the sequel we restrict ourselves to the 
explicit calculation of is, 7, a and the critical point. We try to obtain these quantities 
through the determination of £_l(A£), m2i,(At)/L 2 and m2L(At)/miL(At) where L 
is the linear size of the 2D lattice and t stands either for p (bond probability) or m 
(cluster multiplicity) . 

The choice of m2L(At)/L 2 deserves a comment. In practice we are interested 
in m2L(At), but this quantity does not converge to any finite value when L goes 
to infinity, whatever t (extensive quantity). We found no rigorous mathematical 
proof for the convergence of the aforementioned quantity to a constant asymptotic 
value with increasing L. It is however trivial to see that m,2L/L 2 = 1 when t = 
p = 0. It is equal to when t = p = 1 (the heaviest cluster is not taken into 
account). Furthermore, numerical tests indicate that this quantity may indeed reach 
an asymptotic value, although this value is reached for larger and larger values of 
L when t (p, m) gets closer and closer to t c . We discuss this point in the sequel. If 
this convergence property is fulfilled, one faces a second point which is the question 
whether the exponent and location of the critical point are the same for m^L and 
iti2l/L 2 when L is infinite. One can convince oneself that the corresponding scaling 
functions Q are not the same. The tests presented below confirm that t c and the 
exponent are in numerical agreement with the known values [17]. Finally, it is easy 
to show that ttiil/L 2 tends to 1 with increasing L. Hence m2Ll^\L reaches a 
constant limit for very large systems and allows to extract in principle the exponent 
(a or a, see below). 

Calculations and the results which will be shown correspond to systems without 
periodic boundary conditions. Comments about calculations made for systems with 
periodic boundary conditions will be made below. 
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3a. Calculations with the bond probability p 

We have determined the correlation length £/,(Ap), the normalized second mo- 
ment m2L(Ap)/L 2 and the ratio rri2L(A.p) /miL(Ap) for L up to 1 000 and different 
values of p in the interval [0.40, 0.48] in steps of Ap = 0.01. For each value of p we 
generated 10 4 events. The critical value is p c = 0.5. Fig. 1 shows the behaviour of 
the observables. As it can be seen, these quantities tend to reach a constant value. 
One may notice that this saturation effect is better realized for smaller values of p 
and more effective for £l than for m,2L- As it can be seen on the figure, it is not ef- 
fectively reached for the highest value of p. The consequences of these observations 
will be seen below. 

We retain the asymptotic values reached by the observables as those which 
correspond to the infinite system and use them in order to get the universal scaling 
functions Q^, Qm 2 L/L 2 an d Qm 2 L/m 1L through (1). Fig. 2 shows examples of 
these quantities obtained for different values of p and represented as functions of 
x = £l/L. Universality requires that all contributions from different p's lie on each 
other and asymptotic saturation requires that the Q's reach a constant value (equal 
to 1) for x close to 0. The universality functions are fit to polynomial in x, see eq. 
(14) in ref. [14]. Finally these quantities are used in order to determine p c , v, 7 and 
a by means of a fit procedure to the quantities relevant to the infinite system, i.e. 



£00 ^ (p c -P) " 

and 

m 2o o/m loo = (Pc-p)~ 1/a 

The quality of the fits is estimated by means of x 2 analyses. Fig. 3 shows 
a typical example of the calculations. In practice we worked with many different 
intervals of p values in the range [0.40, 0.48] over which we averaged. We retained 
those results which correspond to the smallest x 2 values, with a confidence level of 
95%. 

Results are shown in Table 1 and Fig. 4. The exponents are calculated from 
systems of different sizes. They are in an overall good agreement with the exact 
values [17]. One notices however that is, 7 and 1/a decreases with increasing L, v 
gets closer to the exact value while 7 and 1/a get away from it. Introducing periodic 
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boundary conditions leads to values which remain very close to those which do not 
take them into account. In practice, the two types of calculations differ essentially 
in the way asymptotic saturation of the observables is reached, the rise being steeper 
when periodic boundary conditions are taken into account. 

3b. Analysis and discussion of the FFS method 

As mentionned above, we have tested the present method by choosing different 
intervals in p lying more or less close to p c . If the interval is taken for, say between 
0.30 and 0.40, the exponents are robust with respect to the value of L from which 
they are extracted. The exponents are however sizably smaller than the exact known 
values [19]. If the interval approaches p c , they come closer to the exact value but 
they are less robust with respect to L, as already mentionned. 

These ascertainments have a common origin. In principle, the universal func- 
tions Qo are obtained from the asymptotic constant limit of O. This limit with 
increasing L is easily reached when p is small. If however p comes close to p c , higher 
values of L are needed in order to reach this limit. It comes out that this limit is 
more easily reached for £l than for m^Ll L 2 and ?712l/toil This reflects on the be- 
haviour of Q^, Qm2L/L 2 an d Qm2L/miL- Indeed, one can observe that for values of 
L up to 1 000, reaches a plateau for the smallest values of x = whereas 
one can hardly see this effect for Q m 2L/L 2 an d Qm2L/miL- This means in practice 
that one needs to determine the considered observables for values of L which are 
sizably larger than L = 1 000. We did not work these cases since it is not our aim 
to work the exponents for 2D percolation which are well known, but rather test a 
method and its efficiency. 

It is possible to find an explanation for the difference in the asymptotic be- 
haviour of the considered observables. It is easy to observe that the rate of con- 
vergence with L decreases with increasing values of the exponents. This fact can 
also be observed in the determination of exponents in ref. [14]. Indeed 7 and \ ja 
are much larger quantities than v. Larger exponents correspond to steeper increase 
of the corresponding observable, hence a larger sensitivity to the precise numerical 
determination of this observable. One may remark that the exponent 7 is sizably 
smaller (~ 1.8) in 3D bond percolation systems. We conjecture that one may need 
smaller values of L in order to reach the asymptotic regime. 
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Sc. Calculations with the multiplicity m 

A confrontation of a percolation model determination of observables like 7712 
and 7712/7711 with experiment [10] is possible if one uses the fragment and particle 
reduced multiplicity m [20] (multiplicity divided by L 2 ). It is numerically easy to 
relate p to an average value of m, < m >. It is then tempting to work out the 
behaviour of the quoted observables £/_,, tt^l/-^ 2 and va2hl^i\L and to confront the 
exponents v, 7 and a with v, 7 and a obtained with the bond probability p. 

We repeated the preceeding calculations with m, following the same procedure 
as described in section 3a. For a given m the events were selected out of a uniform 
distribution in p. For each value of m we generated 2 500 events. The same remarks 
as those presented above are valid in the present case. However, as it can be 
seen in Table 1, the corresponding exponents v, 7 and a are no longer numerically 
compatible with u, 7 and a. In fact, v and 7 are systematically smaller than v 
and 7 whereas a is larger than a. This result contradicts the findings presented in 
ref. [12]. One can convince oneself that the discrepancy cannot be attributed to 
the convergence trouble quoted above, this discrepancy remains for v corresponding 
to £l, for which the saturation problem is much less acute than for 7 and a. In 
fact, this result is not really surprising for us. Indeed, the relation between p and 
< m > is definitely not linear in the neighbourhood of p c as it seems to be for 3D 
percolation on Fig. 7 of [12]. 

The calculations have been repeated for systems with periodic boundary con- 
ditions. As before, the values of the exponents lie within the error bars obtained in 
the case of open boundaries. 

4. Summary, discussion and conclusions 

We have used a rather simple method relying on the finite size scaling (FSS) 
assumption which has earlier been tested on the 2D and 3D Ising model [14]. The 
method allows the rigorous determination of critical exponents characterizing second 
order phase transitions in genuine infinite systems. 

The method is numerically simple, in principle robust, and should deliver expo- 
nents with a very good precision. We verified that it is indeed possible to get values 
which lie within a few percent of the exact value with a rather modest amount of 
work. The central point lies in the fact that one must determine the asymptotic 
value of the corresponding observables as cleanly as possible. This means that it is 
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necessary to determine these observables for larger and larger systems. The maxi- 
mum size corresponding to the asymptotic (infinite size) regime depends on the 
observable itself and empirical experience shows that it depends upon the magni- 
tude of the exponent itself. There exists of course methods which give much more 
precise exponents [17,18] but they are numerically much more sophisticated than 
the present one. 

In the present work we have concentrated on the academic example of 2D 
bond percolation and seen that the exponents depend on the relevant parameter 
(bond probability p or cluster multiplicity m) which is used. One may ask oneself 
whether and how the FSS method could be applied to a realistic physical case such 
as nuclear fragmentation. As we already said in the introduction, the determination 
of critical exponents requires rigorous information about the infinite system. This 
shows that the determination of exponents for nuclear systems may be quite difficult. 
Indeed nuclei are very small (number of coexisting nuclei ~ 300 — 400 in heavy ion 
collisions), unlike systems found in condensed matter physics where surface effect 
can be made negligible if the considered samples are large enough. 

The necessary reference to the infinite system also shows that there exists 
no direct way to extract exponents from experimental data. One way one may 
think about consists of working out a model which reproduces the experimental 
observables obtained for fixed finite sizes L. One may then extrapolate the model 
to L = oo in the FSS framework and determine the exponents the way which was 
described above. If the observables can be determined experimentally for several 
sizes (f.i. Li and L 2 ) and if they are in agreement with the model then, using the 
notations of section 2, these observables must verify the scaling relation 

O e ^(At)/OT»(At)=Qo(x {L u t))/Q {x{L 2 ,t)) 

If this is the case one may conclude that the system is critical in the infinite limit 
and the exponents determined in this limit effectively characterize the criticality of 
the system. 

Whether such tests can be performed or not in practice is an open question 
which is strongly related to the precision with which relevant observables can be 
experimentally determined. 
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Table caption 

Table 1 : Location of the critical point and exponents calculated for different 
linear sizes L of the 2D percolation system using the FSS assump- 
tion (1). The numbers in parentheses correspond to the estimated 
error on the last figure of the given value. 
The exact values are : 

p c = 0.5 ; 7 = 2.389 ; v = 1.333 ; a = 0.396 



L 


Pc 


7 


V 


a 




40 


0.505(5) 


2.43(6) 


1.42(5) 


0.40(2) 




80 


0.503(5) 


2.43(4) 


1.40(6) 


0.40(2) 


a) 


120 


0.501(3) 


2.34(7) 


1.35(3) 


0.42(3) 




40 


0.505(5) 


2.41(8) 


1.38(7) 


0.41(4) 




80 


0.500(2) 


2.41(5) 


1.39(4) 


0.43(3) 


b) 


120 


0.500(2) 


2.28(8) 


1.35(4) 


0.44(3) 




L 


m c 


7 


V 


a 




40 


0.088(5) 


2.08(4) 


1.26(7) 


0.46(2) 




80 


0.097(3) 


2.01(7) 


1.17(7) 


0.49(3) 


c) 


120 


0.098(3) 


1.98(4) 


1.17(3) 


0.50(2) 




40 


0.100(7) 


1.93(7) 


1.15(7) 


0.51(5) 




80 


0.097(3) 


2.00(6) 


1.14(5) 


0.47(3) 


d) 


120 


0.102(3) 


1.97(5) 


1.09(5) 


0.51(3) 





bond probability p with open boundaries 

bond probability p with periodic boundaries 

reduced multiplicities m with open boundaries 
reduced multiplicities m with periodic boundaries 
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Figure caption 



Fig. 1 : Evolution of £ and 1712/ L 2 as a function of the linear size L for different 
values of p. 

Fig. 2 : Behaviour of the universal functions Q related to £ and 1712 / L 2 as a function 
of x for different values of p. 

Fig. 3 : Typical fit of the second moment represented as a function of (p c — p) 
on logarithmic scale for L = 40. The dispersion generated through the 
simulations is smaller than the size of the points. 

Fig. 4 : Values of the critical exponents from bond probability calculations ex- 
tracted from the simulations for different values of L. The exact known 
values are indicated by the lines. 
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